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Abstract. 

Q . The origin of the lowest mass free-floating objects - brown dwarfs and planetary- 

?— i ■ mass objects - is one of the major unsolved problems in star formation. Establishing 

c/3 \ a census of young substellar objects is a fundamental prerequisite for distinguishing 

d . between competing theoretical scenarios. Such a census allows us to probe the initial 

mass function (IMF), binary statistics, and properties of accretion disks. Our SONYC 
(Substellar Objects in Nearby Young Clusters) survey relies on extremely deep wide- 
field optical and near-infrared imaging, with follow-up spectroscopy, in combination 
■ with Spitzer photometry to probe the bottom end of the IMF to unprecedented levels. 

Here we present SONYC results for three different regions: NGC 1333, p Ophiuchus 
and Chamaeleon-I. In NGC 1333, we find evidence for a possible cutoff in the mass 
' function at 10-20 Jupiter masses. In p Oph we report a new brown dwarf with a mass 

\ close to the deuterium-burning limit. 

o 

> : 

•'■h . 1. Introduction 

The origin of the stellar initial mass function (IMF) is one of the major issues in as- 
trophysics. The low-mass end of the IMF, in particular, has been subject of numerou s 



observational and theoretical studies over the past decade (see iBonnell et al.l 120071) . 
Star forming regions and young clusters harbor a large population of free-floating ob- 
jects below the stellar-mass boundary, found to be at ~0.075 M Q . These objects include 
brown dwarfs (BDs), but also a population of objects with masses comparable to those 
of massive planets (we refer to these objects as "planetary-mass objects" or PMOs). De- 
spite the large advances in our understanding of substellar objects over the past decade, 
some of the most i mportant questions still remain unanswered. The origin of BDs and 
PMOs is not clear dWhitworth & Goodwin! 120 05): competing scenarios include t urbu 



lent f ragmentation (IPadoan & Nordl und 2004), ejection fro m multiple systems (Bate 



2009), and ejection from fragmenting protoplanetary disks dStamatellos & Whitworth 



2008). The shape of the IMF at very low masses is the subject of an ongoing debate in 
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the literature ( Bonnell et al.ll2007l : IChabrieril2003l) . In nearby star forming regions the 
total number BDs r elative to the number o f low-mass stars varies between 3 and 8 with 
large uncertainties (lAndersen et al.l l2008). The number of PMOs and its dependence 
on environment is eve n more uncertain. The IMF c ould be still risi ng below 0.015M o 
(iCaballero et al]|2007r) . or declining in this regime (ILucas et all 120051) . A cutoff in the 
mass function has not been observed yet. 

SONYC - Substellar Objects in Nearby Young Clusters - is an ongoing project 
to provide a complete census of the brown dwarf and planetary mass object population 
in nearby young clusters, and to establish the frequency of substellar mass objects as 
a function of cluster environment. The resulting catalog of substellar mass candidates 
will provide the basis for detailed characterization of their physical properties (disks, 
binarity, atmospheres, accretion, activity). The primary means of identifying candi- 
dates is broad-band imaging in the optical and the infrared, thus aiming to detect the 
photosphere. The survey is also combined with the 2MASS and Spitzer photometry 
catalogs. Photometric selection results in large samples of candidates and requires ex- 
tensive spectroscopic follow-up to asses the real nature of the objects. Our observations 
are designed to reach limiting masses of ~ 0.005 M Q , well below the deuterium-burning 
limit at 0.015 M Q , and thus require us of 4- to 8-m-class telescopes. By probing several 
star forming regions we want to probe for environmental differences in the frequency 
and properties of substellar objects. 

In this contribution, we summarize the results delivered in the framework of SONYC 
over the past two years. We hav e surveyed three star forming regions: NGC 1333 
(IScholz et al.ll2009l) . p Ophiuchus dGeers et al.ll2010l) . and Chamaeleon-I (Muzic et al., 
submitted to ApJ). 



2. NGC 1333 



NGC 1333 is a clust er in the Perseus star f orming complex, w ith an age of ~1 Myr, a 
distance of ~300pc dde Zeeuw et al.lll999l : lBelikov et al.l l2002). and moderate extinc- 
tion. The spatial coverage of our survey in NGC 1333 is 0.25 deg 2 (see Fig[TJ). 



2.1. Photometry and Spectroscopy 



Imaging of NGC 1333 was per formed at the Subaru telescope. We used the Suprime- 
Cam wide-field optical i mager ([ Miyazaki et al. 2002) to obtain im ages in SDSS filters 
i' and z', and MOIRCS (Suzuk i et al.ll2008l : llchikawa et al.ll2006l) for near-infrared J 



and K$ observations. Our data are complete down to /'=24.7, z'=23.8, /=20.8, and 
^5 = 18.0 mag. In terms of object masses for members of NGC 1333, this roughly 
corresponds to mass limits of O.O OSMp, for Ay <10 and 0.004MQ for Ay <5, based on 
the COND03 (iBaraffe et al.ll2003l) and DUSTY00 (IChabrier et al.ll2000h evolutionary 
tracks. We used MOIRCS again to carry out multi-object spectroscopy for 53 sources 
in NGC 1333, selected on the basis of their broad-band colors that are in the range as 
expected for substellar cluster members. This sample does not show any bias in spatial 
coverage or optical/NIR colors with respect to the full photometric candidate sample. 
The wavelength coverage for MOIRCS low-resolution spectroscopy includes H- and 
^-bands. 
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Figure 1. Spatial coverage of our survey in NGC 1333. Crosses are all candidates 
selected from optical photometry with i' >18.5, diamonds the spectroscopica l ly con - 
firmed substellar members, triangles the confirmed BDs from Wilk ing et al.l ([2004). 
The plusse s are YSO cand i dates from Spitzer data without spectroscopic confirma- 
tion from lGutermuth et al.l (120081). The underlying co ntours show the distribution of 
13 CO from the 7=1-0 map by iRidge etaT] d2006l) . The dotted and dashed lines 
show the area cov ered in our optical and NIR survey, respectively. Image from 
IScholz et al.l d2009t) ; reproduced by permission of the AAS. 



2.2. Results 

Our goal is to identify young sources wi th effective temperatures at or below 3000 K. 
As outlined in detail in lScholz et al.1 (120091) . these objects show a characteristic spectral 
shape in the near-in frared. In particular, their spectra have a clear peak in the H-band 
(ICushing et alJl2005h . This feature is caused by water absorption on both sides of the 



H-band. The depth of this feature depends strongly on effective temperature. While the 
H-band peak appears round in old field dwarfs, it is triangular in young, low-gravity 
sources. In addition, young brown dwarfs are expected to have flat or increasing /ri- 
band spectra with CO absorption bands at A > 2.3 /urn. To estimate the effective tem- 
peratu res of the BP cand idates we perform fitting of model spectra from the DUSTY 
series ( Allard et alj|200lh to our observed spectra. We confirm 19 sources with effective 



temperatures of 3000 K or lower. Combined with the clear indications for youth, we 
classify them as probable substellar members of NGC 1333. 

Combining our survey results with previous studies, the current census of spectro- 
scopically confirmed BDs in NGC 1333 is 33. The ratio of substellar to stellar members 
with masses below IM Q is 1.5 ± 0.3, lower by a factor of 2 - 5 than in all other previ- 
ously surveyed regions. Thus, NGC 1333 clearly shows an overabundance of BDs. On 
the other hand, the cluster shows lack of PMOs. The low-mass limit of the confirmed 



4 K. Muzic, A. Scholz, V. C. Geers et al. 




Figure 2. Mass distribution of BDs in NGC 1333 (hatched histogram) and the 
deficit of planetary-mass objects. The three labeled data points show the predicted 
number of PMOs in NGC 1333 based on the surv eys in cr Ori , the O NC, and in Cha- 
I. The error bars correspond to lcr. Image from lScholz et a l. (2009); reproduced by 
permission of the AAS. 



BDs is 0.012 - 0.02M o , but the completeness limits are at significantly lower masses. 
Scaling from literature results in other regions (cr Ori, ONC, Cha-I), we would expect 
to find 8-10 PMOs, but we find none (see Fig [2]). This indicates a cutoff in the mass 
spectrum around the deuterium-burning limit in NGC 1333. 



3. p Ophiuchus 



p Oph is one of the closest (d = 125 + 25 pc. Ide Geus et al.llT989h regions of active star 
formation, p Oph cluster is not as compact as the first SONYC target NGC 1333, and 
exhibits extremely high and variable levels of extinction. The main cloud, L1688, is 
a den se molecular core, with visual extinction up to 50 - 100 mag ( Wilking & Ladal 



1983), hosting an embedded infrared cluster of around 200 stars, inferred to have a 



median age of 0.3 My r, and surrounded b y multiple clusters of young stars with a 
median age of 2. 1 Myr (|Wifking et al. 2005 and references therein). 



3.1. Photometry and spectroscopy 

For p Oph we used the same instrumental setup as for NGC 1333 (Section 12.11) . to 
obtain broad-band images in i' , J and K$, and HK spectroscopy. Completeness limits 
are placed at i'=24.2, 7=20.6, and #s = 17.8. This corresponds to mass limits of 0.004 
- 0.1M o in the /'-band, and 0.001 - 0.007 M Q in the /-band, for extinction of Ay - 5 - 




Figure 3. Spatial distribution of sources in p Oph. Contours are Ay = 5, 10, 
15, 20, 25, 30, as derived from 2MASS by the COMPLETE project. Dotted line: 
/-band imaging coverage; dashed line: J and /Tj-band imaging coverage; dash-dot 
line: MOIRCS spectroscopy fields. iJK$ catalog sources indicated with +; UK$ 
catalog selected BD candidates with follow-up MOIRCS spectroscopy indicated 
with x; spectroscopically confirmed low mass YSOs are indicated with red large 
squares. Candidate BDs selected f rom JKs + Spitzer p hotometry indicated with cir- 
cles. Previously k nown BDs from Wilking et al. ( 2008) and class II / III sources from 
iNatta et alj d2006h are indicate d with small blue squares and diamonds respectively. 
Figure from lGeers et al.l d2010h ; reproduced by permission of the AAS. 



15 (based on the COND03 and DUSTY00 evolutionary tracks). Photometric catalogs 
were cross-correlated with the existing Spitzer data. The spatial coverage of our i'+JK$ 
survey is 0.171 deg 2 (see Fig [3]). 

From the optical and near-infrared photometry, 309 objects were selected as can- 
didate substellar cluster members. 58 of these objects, and 1 additional previously 
known BD candidate, were targeted for follow-up spectroscopy. Based on multi-object 
spectroscopy, using the water absorption features in the H-band, 1 of the 58 new can- 
didates was confirmed as a substellar mass object with T e s = 2500 K (see left panel 
in Fig 01). From MOIRCS, 2MASS, and Spitzer photometry a sample of 27 sources 
with mid-infrared color excess and near-infrared colors indicative for substellar mass 
sources with disks are identified. Of these, 1 1 are previously spectroscopically con- 
firmed brown dwarfs, while 16 are newly identified candidates. Based on present day 
surveys of the stellar and brown dwarf populations, the ratio of substellar to stellar 
sources in p Oph is derived to have an upper limit of 5 - 7, in line with other nearby 
young star forming regions. Census of substellar objects in p Oph, based on current 
existing surveys, is likely highly incomplete, due to the variable and high extinction. 
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Figure 4. Left panel: MOIRCS spectra (solid) of candidate substellar objects in 
p Oph, with best fit reddened model (dashed) with an T e g of 2500 K, 3100 K, and 
3400 K; offsets were applied for clarity. Figure from lGeers et alj d20 1 Ot) : reproduced 
by permission of the AAS. Right panel: Examples of the model fitting for some of 
the previously known Cha-I candidates. In black are shown t he spectra obtain ed by 
VIMOS/VLT, with the best-fit models in blue (AMES-Dustv: lAllard et al1l200lh . 



4. Chamaeleon I 



At a distance of ~ 160pc (Whit tet et al.lll997l). Ch amaeleon I (hereafter Cha-I) harbors 



a rich population of known YSOs (ILuhm an 2007). Compared to the other two SONYC 



targets, Cha-I is less compact than NGC1333 and has less extinction than p Oph. The 
spatial coverage of our survey in Cha-I is 0.25 deg 2 . 

4.1. Photometry and spectroscopy 

This study is primarily based on data from the European Southern Observatory, ob- 
tained in runs 078.C-0049 and 382.C-0174. The fi rst campaign provided deep opti- 
cal iz imaging (VIMOS /VLT; Ee Fevre et all 120031) . and near-infrared JK$ imaging 



(SOFI/NTT: iMoorwood et al.l ll998). In our imaging survey, we reach completeness 
limits of 23.0 in 7, 18.3 in J, and 16.7 in K$. This corresponds to mass limits of 
(0.002 - 0.01)M o , for A v - - 10, based on DUSTY00 and COND03 models. From 
the optical dataset we selected a list of candidate members which were observed us- 
ing multi-object spectrograph VIMOS/VLT, using the low resolution red grism (5500 
- 9000 A). Additionally, we made use of the archival Spitzer images for our target 
regions. 

4.2. Results 



From the optical photometry, 142 objects were selected as candidate substellar cluster 
members. 60 of these objects have been observed in the spectroscopy campaign, how- 
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ever, due to the observed limit for spectroscopy at I«21, only 18 of these candidates 
have spectra that can be used for classification. In addition, we obtained spectra for 
more than 200 randomly selected objects in the VIMO S field-of-view. Red portion o f 



M-dwarf spectra is dominated by molecular features (IKirkpatrick et ail 1199 U 11995). 



which allows relatively simple preliminary selection of candidates based on visual in- 
spection. To determine the effective temperature of the re maining candidate objects we 
performed spectral fitting using the AMES-dusty models (lAllard et aill200lh . We iden- 
tify 13 objects consistent with the spectral type M, of which 11 are previously known 
M-dwarfs with confirmed membership in the cluster (see right panel in Fig|4]). The two 
newly reported objects have effective temperatures consistent with masses above the 
substellar limit. 

Based on the results of our survey and combined with the numbers of substellar 
objects from the literature, we estimate that the number of the missing low-mass mem- 
bers down to ~ 0.008 M Q for Ay < 5 in Cha-I is < 7, i.e. < 3% of the total number of 
members according to the current census. We might, however, still miss objects with 
lower masses, and objects at higher extinctions. 



5. Conclusions 

It is clear that the census of BDs and PMOs in most star forming regions is still incom- 
plete. Based on the existing data we can conclude that: (a) there are hints of regional 
differences in the mass function at the very-low-mass end, and (b) only a combination 
of different search techniques can provide a robust picture of the substellar population. 
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